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ABSORPTION  AND  TRANSFORMATION  OF 
ELECTROMAGNETIC  WAVES  IN  THE  VICINITY  OF 
ELECTRON  CYCLOTRON  HARMONICS 

I.  INTRODUCTION 

Recent  advances  in  the  development  of  multi-kilowatt  RF  sources  in  the  10  10  Hz  fre- 
quency range1  have  led  to  interest  in  the  possibility  of  electron  cyclotron  heating  in  tokamaks. 
Electron  cyclotron  frequency  heating  in  plasmas  has  been  demonstrated  in  smaller  (lower  mag- 
netic field)  devices2  and  preliminary  experiments  in  tokamaks  have  yielded  encouraging 
results.3  In  this  paper  we  will  investigate  the  linear  propagation  and  absorption  of  electromag- 
netic waves  in  the  vicinity  of  cyclotron  harmonics,  that  is  regions  of  space  where  the  wave  fre- 
quency is  an  integer  multiple  of  the  electron  cyclotron  frequency. 

Previous  analyses  have  derived  electron  heating  rates  by  integrating  the  equations  of  mo- 
tion for  the  electrons  in  given  electromagnetic  fields.  This  method  neglects  the  effects  of  the 
perturbed  charge  and  current  densities  in  determining  the  fields  themselves  and  can  lead  to  er- 
ror. 

Our  approach  is  to  analyze  the  propagation  of  electromagnetic  waves  in  an  inhomogene- 
ous magnetic  field.  We  assume  the  scale  length  for  variation  of  the  magnetic  field  and  electron 
density  is  greater  than  both  the  electron  gyroradius  and  the  vacuum  wavelength.  Thus,  a 
WKB  analysis  is  expected  to  be  accurate  except  in  the  immediate  vicinity  of  a cyclotron  reso- 
nance. Near  the  resonant  surface  a full  wave  theory  is  required  that  will  connect  the  WKB 
Manuscript  submitted:  August  5,  1977. 
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solulions  valid  on  either  side  of  the  resonance.  The  full  wave  solution  yields  transmission, 
reflection,  absorption  and  linear  mode  conversion  coefficients  from  which  the  electron  heating 
rate  can  be  determined.  Such  a theory  is  presented  in  this  paper. 

The  organization  of  the  paper  is  as  follows.  In  Section  II  we  briefly  review  the  propaga- 
tion of  electromagnetic  waves  in  a cold  plasma.  In  Section  III  we  present  a simplified  analysis 
of  wave  absorption  at  a cyclotron  harmonic  based  on  WKB  theory.  That  is  we  calculate  wave 
absorption  by  integrating  the  imaginary  part  of  the  wave  vector  across  the  resonance.  Strictly 
speaking,  this  method  is  invalid  if  dkx/dx  > kx  where  kx  is  the  wave  number  and  x the  dis- 
tance from  the  resonance.  However,  this  procedure  is  found  in  later  sections  to  give  the  correct 
transmission  coefficient. 

In  Section  IV  we  derive  from  the  Vlasov  equation  an  expression  for  the  perturbed 
current  density  resulting  from  arbitrary  wave  fields  near  a cyclotron  resonance.  In  Section  V 
and  Section  VI  we  use  the  result  of  Section  IV  to  analyze  the  propagation  of  the  extraordinary 
wave  near  the  second  harmonic  and  the  ordinary  wave  near  the  first  harmonic  respectively. 
Section  VII  contains  our  conclusions  concerning  the  applicability  of  these  mechanisms  to  the 
heating  of  tokamaks. 

II.  REVIEW  OF  COLD  PLASMA  WAVE  PROPAGATION 

There  are  two  electromagnetic  waves  that  propagate  in  a cold,  magnetized  plasma.  If  these 
waves  propagate  in  a direction  nearly  perpendicular  to  the  applied  magnetic  field  they  are  called 
the  ordinary  and  extraordinary  waves. 

The  ordinary  wave  is  characterized  by  the  alignment  of  the  wave  electric  field  nearly 
parallel  to  the  applied  magnetic  field.  For  purely  perpendicular  propagation  the  dispersion  rela- 
tion for  the  ordinary  wave  in  a cold  plasma  (immobile  ions)  is  given  by. 
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(1) 

where  kx  is  the  wave  vector  perpendicular  to  B0,  the  applied  magnetic  field,  c is  the  speed  of 
light,  oj  is  the  wave  frequency  and  w2  — 4ire2n0/me  is  the  square  of  the  plasma  frequency. 
Notice  that  Eq.  (1)  does  not  involve  B0  explicitly.  This  is  a consequence  of  the  parallel  align- 
ment of  the  wave  electric  field  and  the  applied  magnetic  field. 


The  extraordinary  wave  is  characterized  by  the  alignment  of  the  wave  electric  field  per- 
pendicular to  the  applied  magnetic  field.  For  purely  perpendicular  propagation  the  dispersion 
relation  for  the  extraordinary  wave  is  given  by, 

k2  c2  (<o 2 - Wp  — 11  2 ) = (o> 2 — w 2 ) 2 — 11  2a> 2, 


(2) 

where  11  = eB0/mec  is  the  electron  cyclotron  frequency.  The  extraordinary  wave  is  resonant 
(i.e.  k±  — °° ) as  w2  — w2  + 11  2,  (this  is  known  as  the  upper  hybrid  resonance)  and  it  is  cut 
off  (kx  — 0)  as  <o 2 •—  o> 2 + 11  io. 


Including  a finite  k.  (parallel  to  B())  modifies  these  dispersion  relations  slightly.  Figure  1 
shows  an  Allis4  diagram  (assuming  k.  is  fixed)  which  details  the  cut-offs  and  resonances  for 
the  two  electromagnetic  modes  in  a cold  plasma. 


We  wish  to  point  out  that  neither  of  the  waves  experience  a resonance  or  cut-off  when 
the  wave  frequency  is  equal  to  the  cyclotron  frequency  and  k.  is  fixed.  Thus  in  an  inhomo- 
geneous plasma  where  the  density  is  nearly  uniform  along  a field  line  cold  plasma  theory 
predicts  that  electromagnetic  waves  will  propagate  through  regions  where  w = 11  (x)  without 
loss  of  energy.  Linear  absorption  of  wave  energy  by  the  particles  is  a kinetic  effect  that  can 
only  be  determined  by  considering  the  behavior  of  the  plasma  on  a microscopic  level  and  it 
must  result  from  the  plasma  having  a nonzero  temperature. 
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III.  WKB  ANALYSIS  OF  WAVE  INTERACTIONS  AT  THE  CYCLOTRON  HARMONICS 


The  dispersion  relations  for  the  ordinary  and  extraordinary  waves  obtained  from  the 
Vlasov-Maxwell  set  of  equations45  indicate  that  these  waves  will  have  resonances  and  cut-offs 
near  each  of  the  cyclotron  harmonics.  If  the  gradients  in  applied  magnetic  field  and  density  are 
gentle  then  a WKB  description  is  expected  to  be  valid  everywhere  in  the  plasma  except  in  the 
immediate  vicinity  of  resonance  or  cut-off. 


In  the  later  sections  of  this  paper  we  will  present  a "full  wave"  calculation  that  is  valid  in 
the  immediate  vicinity  of  a cyclotron  harmonic  resonance.  The  results  of  our  analysis  is  a solu- 
tion for  the  wave  electromagnetic  field,  valid  in  the  resonance  region,  which  asymptotically 
matches  on  to  the  WKB  solutions  far  from  the  resonance.  Hence,  we  obtain  a set  of  connec- 
tion formulas  which  determine  the  fractions  of  wave  energy  that  are  transmitted,  reflected,  ab- 
sorbed, and  converted  to  other  wave  types. 


It  is  instructive  to  apply  WKB  theory  first  to  the  resonance  region  to  determine  whether 
a wave  will  suffer  significant  energy  loss  in  propagating  through  the  resonance.  WKB  theory  of 
course  is  not  completely  accurate  in  that  it  will  not  predict  reflection  or  mode  conversion 
although  it  is  surprisingly  accurate  in  determing  the  wave  energy  transmission  coefficient.6 


A.  Ordinary  Wave 


The  dispersion  relation  for  the  ordinary  wave  propagating  perpendicularly  to  an  applied 
magnetic  field,  obtained  from  the  Vlasov  Maxwell  set  of  equations,  is  given  by, 

iifi  _ , - r , Jn  ( )v.  6/ 

to2  „ to  J v to  - nil  av.  ’ ,,, 


where  /(v ) is  the  velocity  distribution  function  for  electrons  and  Jn  is  the  ordinary  Bessel 
function. 
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To  determine  the  total  amount  of  wave  energy  damping  in  the  resonant  region  we  as- 
sume that  11  is  a weak  function  of  position,  11  -=  11  (x)  and  integrate  the  imaginary  part  of 
kL(x)  as  determined  by  Eq.  (3)  through  the  resonant  region.  The  effect  of  the  resonance  will 
be  treated  as  a perturbation  to  the  cold  dispersion  relation.  For  example  near  the  first  cyclotron 
harmonic  we  obtain  from  Eq.  (3)  by  expanding  the  Bessel  functions  for  small  arguments 

k2  c2  2 k2\2 

-1  = J _ 'IJL.  , + _i_^£ “ . 

oj  2 cu 2 4(1  2 to  — 11  (4 ) 


Taking  11  (x)  to  be  11  (x)  = wjl  — — J where  L is  the  scale  length  for  variation  of  the  mag- 
netic field,  we  find 


2kkx  (x)c2 

- Im 

°>2r 

1 /c2v2 

1 Klrwe 

i 

» — im 

OJ  2 

OJ  2 

4 cu 2 

X 

L 

+ /€ 

where  kxr  (kit)  is  the  real  (imaginary)  part  of  kx,  and  e is  a positive  infinitessimal  resulting 
from  the  causality  condition  Im  (to)  > 0.  The  wave  energy  transmission  coefficient  7"is  given 
by 


r=expt  I dx2kXl  (x) | =exp(-2rrrj  ) 


where 


2 2)1/2  t 

L <u/>  . wp  \ w Te 

Va  " T ~T  1 21  T 2 

4 iu  cu  ) c mfci 


and  Te  = - The  wave  transmission  coefficient  at  first  glance  would  appear  to  be  nearly  unity  (no 

wave  absorption)  owing  to  the  smallness  of  Te/(m  c2).  However,  the  quantity  Lu>/c  can  be 


quite  large  indicating  small  wave  transmission  and  potentially  high  wave  absorption.  In  a plas- 
ma with  the  following  parameters  L - 130  cm,  B = 35  KG,  Te  - 2 KeV, 
n =0.5  x 10,4cm  ~3,  we  find  T — 5.3  x 10  ~3  indicating  negligible  wave  transmission. 

If  the  same  analysis  is  applied  to  the  ordinary  wave  at  the  second  harmonic  we  find  a 
similar  expression 

T = exp  ( -2irt)b), 

where 


Vb 


Lai 

c 


1 - 


3/2  2 

WP 


(6) 


the  expression  for  in  Eq.  (6)  is  smaller  than  that  in  Eq.  (5)  for  na  by  a factor  of  order 
Te/mec 2 thus  the  transmission  coefficient  at  the  second  harmonic  is  nearly  unity  and  wave  ab- 
sorption and  transformation  are  not  expected  to  be  important  in  this  case. 

B.  Extra  Ordinary  Mode 

For  purely  perpendicular  propagation  the  dispersion  equation  for  the  extraordinary  mode 
is 


\_kW_ 

T-  ^ 


~ K. 


k}  c 2 
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1 "X 
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klc 2 


Y ~ Kn 


klc 2 


~ K, 


E _ 
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where  E±  =2  ~]/2(Ex  ± iEv ),  the  wave  electric  field  has  Cartesian  components 
E = (Ex.  Ey.  0),  the  wave  vector  has  components  k = (kx,  0,  0)  and  the  applied  magnetic  field 
is  in  the  f direction.  The  matrix  elements  Kr,  and  Kri  contain  the  response  of  the  electron 
plasma  and  are  given  by,4 
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where  b = kxve/il,  (Tp  = y wv2),  and  we  have  assumed  a Maxwellian  distribution  of  elec- 
tron velocities.  For  our  purposes  here  it  is  only  necessary  to  have  approximations  for  Krl  and 
Kr/  valid  for  A:2v2  < < <u2. 

These  approximations  are  found  to  be,4 

*,/  = 1 - 


a2 

kxC2  5 

If 

3 

k*c* 

1 ± p 

a,2  1 ± 2/3 

32 

n4 

1 ± fi 


Kr,  = 


]_  2 
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kxc 


1 -p‘ 


(9a) 


(9b) 


where  a2  = a>2/a>2,  fi  = ll/to,  and  fi  = 2 Te/{mec2).  The  first  terms  of  Eq’s.  (9a)  and  (9b) 
are  correct  to  first  order  in  fi.  When  the  dispersion  equation,  Eq.  (7),  is  solved  near  the  first 
harmonic  resonance  due  to  cancellations  it  is  necessary  to  retain  terms  that  are  second  order  in 
fi.  Thus  we  have  included  the  last  term  in  Eq.  (9a)  which  is  second  order  in  fi. 

The  analyses  of  the  preceeding  section  can  be  carried  out  in  a similar  fashion  to  deter- 
mine the  wave  energy  transmission  coefficient  for  the  extraordinary  wave  at  the  first  and 
second  electron  cyclotron  harmonic.  We  find  T = exp  ( —2 irn)  where  at  the  first  harmonic 

I 3/2 


and  at  the  second  harmonic. 
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where  a2  = <o2/oj2. 

In  contrast  to  the  case  of  the  ordinary  wave  the  wave  transmission  for  perpendicular  pro- 
pagation is  nearly  unity  at  the  first  harmonic  and  possibly  low  at  the  second  harmonic. 

The  WKB  treatment  of  wave  absorption  at  a cyclotron  harmonic  for  arbitrary  angle  of  in- 
cidence, i.e.  arbitrary  k:,  can  be  formulated  in  the  following  way.  We  wish  to  solve  the  matrix 
equation  D (k,  x)  • E =0  for  kx  (x ),  where  D (k,  x)  is  the  local  dielectric  tensor  D = 
kvkr  — (tu2/c2)  1 + (47t/(u/c2  )<r  and  a is  the  conductivity  tensor.4  In  the  limit 
v2/c2  < < 1,  D may  be  approximated  by  the  cold  dielectric  tensor  D0  plus  a resonant  term 
D|,  which  results  from  finite  Larmor  radius  effects.  The  tensor  D ^ will  only  be  important  near 
a resonance  (it  contains  terms  that  vary  as  x -1  where  x = 0 is  the  location  of  the  resonance) 
and  we  will  treat  /), , as  a perturbation.  The  solution  for  the  perpendicular  wave  number 
kx(x),  and  the  electric  field,  E(x),  will  differ  slightly  from  their  values  predicted  by  cold 
theory  kx  = k0  + kj , E = E0  + E, , where  D0  (k0,  x)  ■ E0  = 0. 


We  Taylor  expand  D in  kx  about  k0  and  multiply  the  dispersion  equation  on  the  left  by 
E0+  and  obtain  an  expression  for  Arj 

*l  <*)  -E  + • D,  (k0.x)  • E0  E0+  .JL-Doiko).  E0 
the  wave  transmission  will  then  be  given  by  T = exp  ( — 2ttt\  ),  where 

oo 

V = — f dxk,  (x). 
rr  i —oo 

Calculations  of  this  type  have  shown  that  the  extraordinary  wave  can  be  absorbed  at  the  first 
harmonic  for  sufficiently  large  kl . 
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IV.  PERTURBED  CURRENT  DENSITY  IN  AN  INHOMOGENEOUS  MAGNETIC  FIELD 

In  this  section  we  begin  the  "full  wave"  analysis  of  electromagnetic  wave  propagation 
near  a cyclotron  harmonic  resonance  that  will  provide  us  with  connection  formulas,  linking  the 
various  WKB  solutions  that  are  valid  far  from  a resonance.  Our  approach  will  be  to  assume  the 
wave  electric  field  can  be  expressed  as  a complex  Fourier  integral,  then  determine  the  response 
of  a plasma  from  the  linearized  Vlasov  equation. 

Substituting  the  given  form  of  the  electric  field  and  the  perturbed  current  density  in 
Maxwell’s  equations  will  yield  an  integral  equation  in  k- space  which  we  will  solve  in  later  sec- 
tions. Once  the  integral  equation  in  k- space  is  solved  for  the  Fourier  representation  of  the  wave 
electric  field  the  connection  formulas  are  obtained  by  asymptotically  expanding  the  Fourier 
inversion  of  the  /c-space  representation  of  the  electric  field. 


Our  model  of  the  plasma  is  the  following.  We  assume  that  the  applied  magnetic  field 
points  in  the  z direction  and  depends  on  x.  The  density  in  general  is  also  a function  of  a- how- 
ever, near  a cyclotron  harmonic  resonance  the  important  spatial  dependence  is  that  of  the  mag- 
netic field  and  hence  the  variation  of  the  density  near  a resonance  will  be  neglected.  We  take 
the  electric  field  to  be  of  the  form, 


dK L 

1 . 2tt 


E(.v)  = J*  E (kx  )exp  (ikxx  + ik7z  — iwt). 


(10) 


where  we  have  neglected  variations  of  the  wave  quantities  in  the  y direction,  and  L is  a con- 
tour in  the  complex  Av-plane  that  will  be  specified  later. 


In  order  to  determine  the  response  of  the  electrons  we  solve  the  Vlasov  equation  using 
the  method  of  characteristics,5  that  is, 
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Li 


where/,  (f0)  is  the  perturbed  (unperturbed)  distribution  function,  q(m)  is  the  particle  charge 


(mass)  and  the  quantities  in  the  integrand  are  evaluated  along  the  unperturbed  orbits  of  the 


charged  particles.  The  applied  magnetic  field  is  assumed  to  have  a linear  dependence  on  x in 


the  immediate  vicinity  of  the  cyclotron  resonance,  B = BA  1 - x/ZJ,  where  L is  the  scale 


length  of  the  magnetic  field  and  isassumedto  be  much  larger  than  either  thegyroradius,  v(,/ii, 


or  the  vacuum  wavelength  c/u>. 


The  particle  orbits  to  lowest  order  in  L 1 are  given  by. 


x'(t)  — x = (vj/O  ) [sin  (fir  -0)  + sin0], 


y'(r)  — y = (vj/ll  ) [cos  ((It  - 0)  - cos  0 ] + - ^ ^ , 


v (r  ) = Vj.  COS  (Or  — 0 ), 


vy  (t ) = Vj.  sin  (fir  — 0 ), 


where  (1  =00[1  — x/ L — Vj./ (fl nL)  sin  0 ] and  11 0 =qB0/(mc).  11  is  the  cyclotron  fre- 


quency of  a particle  with  its  guiding  center  located  at  x + (vx/!l  0 ) sin  H.  It  may  appear  at  first 


that  it  is  permissible  to  neglect  the  (vj/1 10L)  sin  H dependence  of  11  owing  to  the  smallness 


of  the  gyroradius  compared  to  the  scale  length.  However,  this  term  is  found  to  have  an  impor- 


tant effect  near  a resonance,  <o  =«110(1  — x/L). 


To  simplify  our  analysis  we  switch  from  fixed  Cartesian  coordinates  to  rotating4  coordi- 


Ex  = 2 “i/2(E.  + E_  ), 


Ey 12  ~V2(E+  —/•_), 


v'  = vx  cos  v ' = Vj.  sin  d>',  d>'  = 0 —Hr, 


I 


where  we  have  neglected  the  diamagnetic  drift  in  Eq.  (lib).  Of  interest  is  J = A-nqmlw, 
where  nv  is  the  perturbed  current  density, 


-,Vjvv>/ 


9v 


Expressed  in  rotating  coordinates  we  have 


(kx,  x)  ■ E (,kx  )exp  [//rxx  + ik.z  — o>r). 


(12) 


where. 


I = - 


2 w 


U 

J \j_d\iduM  J dr  exp  [/i(sin  (Or  - 0)  +sinfl)  -wr  + ik,ur] 


df0 

v,  — — exp  (/(It) 

ovj. 

^ j' 

Vi.  — — exp  (/II  r —2 iO) 


dvj. 


bf0 

Vi  — — exp  [/( 20  — 11  t )J 

dVi 

9/0 

vx- — exp  ( -/(It) 

9 Vi 


2 1/2u exp  (/(ll r — «)]  2 l/2u  exp  ( —/((It  -0)1 


■ "i  1/2, 


fi/o 


dvi 


9vi 


2,/2Vi4^exp(,«) 

9m 


21/2Vi^Lexp(-/«) 

9m 


2m 


% 
9 m 


and  b = kx\xH\.  Here,  the  x dependence  of  £ is  a result  of  the  x dependence  of  11,  otherwise 
the  conductivity  tensor  I is  the  same  as  for  a homogeneous  plasma.4 


We  adopt  the  notation  of  Allis4  for  the  various  elements  of  the  tensor  £ 


I = 


£ 


~n> 


pr 


*1 

£ 


pi 


We  now  evaluate  I,, 

r 9/„  /•  » 

Zr  “ — -r— J ViiA'if/WrfMVi  — — J r/r exp  { / (/> (sin  (Hr  -0)  + sin0)  + (11  + A.m  — w)t]1 


9vi 


11 


Using  the  Bessel  identity  exp  (ib  sin  (Hr  — 0))  = T.nJ„  (b)  exp  (/>»  ( 12  r - U))  the  integra- 
tion over  t may  be  performed  to  yield, 


I, 


r , 9/0  J„  (i)exp  [/ (6 sin  H - nU)] 

\ v}dv±dt)du-~l„  

J 1 0Vi  " (»  + 1 )ll  +k.u  - w 


(13) 


Here  12  is  a weak  function  of  x,  v,.,  and  d.  We  will  treat  12  as  a constant  everywhere  except 
when  it  appears  in  a resonant  denominator.  That  is,  if  a,  = (m  + 1)120  we  will  retain  the 
x,  vA,  and  H dependences  of  12  in  the  denominator  of  the  m,h  term  of  the  sum  in  Eq.  (13),  and 
everywhere  else  we  will  replace  12  by  12 0.  With  this  approximation  divides  naturally  into 
two  parts;  a part  which  depends  on  x,  Ifl , (the  m'h  term  of  the  sum)  and  a part  which  is  in- 
dependent of  x,  i.ro  (the  remaining  terms  of  the  sum).  ir0  may  be  evaluted  in  the  same 
fashion  as  one  would  evaluate  an  element  of  the  homogeneous  plasma  conductivity  tensor.  The 
evaluation  of  proceeds  as  follows: 

■2  - d/0  Jm  (i)Lexp  [/' (b  sin  H -mb)} 

dvA  x-  + (vx/12 0 ) sin  b + Lk.u/ui  + ieL' 


= + 


2a,2 


J v 2 d\x  dudb 


where  t is  a positive  infinitessimal  resulting  from  the  causality  condition. 

The  a dependence  of  f rl  can  be  expressed  as  the  inversion  of  a Fourier  transform, 

oo 

tr i = f dkt  exp  (/7c,  x)Zrl  (kx,  A,  ) 


where 


(**.  *1  ) “ TT  f 'l^dud*  ^2-  LJm  (6) 


2,. 2 


Vo 


x exp  I + ik]  ( (va/12 0 ) sin  b + Lk.u/w)  + i(b  sin  « — wfl)]. 
Using  the  Bessel  identity  we  find, 

TT  Ixitl  i r , 9/„ 

= - — J \}d\i.du  - — exp  (/A,  Lk.u/w  )Jm  ( b)Jm  (b  4-  6,  ) 

U)  1 


•a&i 


12 


! 

t 

H 


where  ft  = Aj  vx /II I,,  can  be  further  reduced  when  f0  is  specified.  If  f0  is  Maxwellian 
(f0  — 77  ~3/2\e~3  exp  ( - v2/v 2 ) ) we  find 

tt<Lo> 2 d 


at2  9m 


m exp 


-p(J,2  + 2ftft,  + 2ft2) 


y M ft(ft  + V 


x exp  ( — a) 


(15) 


1 


where  ft  = kxve/tt  0,  b\  = klve/il0,  and  a = y (A]  k:L\ze/u> ) 2.  If  /0  has  a Maxwellian 

dependence  on  perpendicular  velocity  and  a Lorentzian  dependence  on  parallel  velocity, 
(/0  — rr  ^v^-1  exp  ( — v2/v2 ) (m2  + v 2 ) _1 ) we  find  the  Irl  is  given  by  Eq.  (15)  with 
a — klk:L\e/<u.  This  latter  form  of  If)  will  be  useful  in  determining  the  effects  of  nonzero  k. 
on  wave  absorption  and  mode  conversion. 

Finally,  it  remains  to  write  down  an  expression  for  the  perturbed  current, 

Jdk * „ 

-r^expOVUr)  H0  (kx ) • E (kx)  +1,  (kx,  x)  ■ E(kx)] 

L 

writing  Ij  in  terms  of  its  Fourier  inversion,  and  changing  variables  of  integration  we  find 


r dkx 

J = J y—  exp  (ikxx) 


2tt 


X 

i0(kx)  E + J*  dk T,  (k',kx  - k')  ■ E(A') 


(16) 


where  the  first  elements  of  and  Ii  have  been  determined  in  this  section.  The  remaining 


elements  can  be  obtained  using  the  same  technique. 


V.  EXTRAORDINARY  WAVE  AT  THE  SECOND  HARMONIC 

The  full  matrix  equation  for  wave  propagation  near  a cyclotron  harmonic  in  an  inhomo- 
geneous magnetic  field  is  obtained  by  substituting  J as  given  by  Eq.  (16)  in  Maxwell’s  equa- 
tions. The  result  is  three  coupled  integral  equations  for  the  three  components  of  E(AV).  To 
analyze  these  equations  we  assume  nearly  perpendicular  propagation  k,e/u>  < < 1.  In  this  limit 
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t 


the  ordinary  and  extraordinary  waves  decouple.  That  is,  the  propagation  of  the  ordinary  wave  is 
governed  by  a single  integral  equation  for  E:(kx)  and  the  propagation  of  the  extraordinary 
wave  is  governed  by  two  coupled  integral  equations  for  E+  (kx)  and  E_  (kx).  The  effects  of 
small  k:  are  retained  in  the  extraordinary  wave  calculation  by  including  k.  in  the  evalution  of 
the  matrix  elements  of  I.  What  we  find  is  that  for  sufficiently  large 
k..\  > > k,c/w  > c/(vekxL)  the  results  of  simple  WKB  theory  become  valid.  This  occurs 
when  the  spatial  width  of  the  cyclotron  resonance  (|oj  — (m  + 1 )fl  (x)|  = k.\e)  is  compar- 
able to  the  perpendicular  wave  length  ( k± ) predicted  by  cold  plasma  theory.  Thus,  departures 
from  simple  WKB  theory  occur  only  for  the  case  of  nearly  perpendicular  propagation. 


We  consider  the  propagation  of  the  extraordinary  wave  near  the  second  harmonic  in  the 
limit  v2  < < c2.  In  this  limit  the  only  resonant  term  is  Irl . The  remaining  nonresonant 
terms  I0  can  be  approximated  by  their  cold  plasma  limits.  Thus  Maxwell’s  equations  reduce  to 


r dkx 

0 = J 

T 27t 


exp  (ik  x) 


1 

*2c2 

2 

o>2 

-1  -!„£+(*,) 


2.2 


1 kxc 


k„k'e2 


2 <o 2 


E - (kx  ) -2iLa2S  f dk'-^—z — exp  [d‘  -d]E+(k') 


r dkx 

0=1  — — exp  (/Ac  jc) 

L 27r 


1 

2 c2 


E+(l<x)  + 


1 

kx  °2 

2 

u>2 

- 1 - I, 


E_ 


(17a) 


(17b) 


where  Iro  = -a2/(l  — /3 ),  I/o  = -a2/(l  + fi),  d = kxk,L\e/w,  and  d'  = k'k:L\e/u>, 
and  we  have  assumed  a Lorentzian  distribution  in  axial  velocity.  Equations  (17a)  and  (17b)  are 
satisfied  when  the  integrands  vanish.  Thus,  we  can  determine  E _(kx)  in  terms  of  E+  ( kx ) 
from  Eq.  (17b)  and  substitute  in  (17a)  to  obtain  the  single  equation, 


(kx)E+(kx)exp(d)/kx  - iL  J dk'G(k')E+ (k')exp  (d')/k'  - 0. 


(18) 


where 
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ttx)  = J k2C2/(o2  - 1 -lr0 


-\k2C2l „2  / \k2c2l»2  -1  -I/o 


G(k')  = 22a  2Sk'2c2/oj 2 

By  differentiating  Eq.  (18)  with  respect  to  kx  a first  order  differential  equation  is  obtained,  its 
solution  is  found  to  be, 


E+  ( kx ) = kxexp  - d + iL  £ dk'G(k')/t  ( k' ) A (kx ) 


which  may  be  inverted  to  yield. 


r krdkx  i / 

E+(x)  =J  — exp  ikxx  — d + iL  \ dk'G{k')/f  (Ar  ) / t(kx). 

L ■»  | / (20) 
Here  we  see  the  quantity  d = Lkxkz\ e/<u  in  the  exponent  of  Eq.  (20)  serves  to  shift  the  x axis 
by  an  amount  ikz\ieL/(D.  So,  we  can  eliminate  d by  making  the  transformation 
x x + ikzxeL/<u.  In  what  follows  we  will  let  d—  0 knowing  that  we  can  recover  the  effects 
of  finite  d by  shifting  the  x axis. 


In  order  to  determine  the  wave  transmission,  reflection,  and  conversion  coefficents  it  is 
necessary  to  asymptotically  evaluate  Eq.  (20)  for  large  positive  and  negative  x.  Before  specify- 
ing the  contour  it  is  helpful  to  see  what  sort  of  asymptotic  behavior  can  occur.  There  are  two 
branch  points  of  £+  ( kx ) both  where  « = 0.  They  are  at 

& i fl2  - 0 ~a2)2 

(O2  p2  + a2  - 1 (2J) 

which  are  the  two  extraordinary  modes  in  a cold  plasma.  In  addition  there  are  two  saddle  point 
contributions  from  saddle  points  at 
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The  steepest  descent  path  makes  an  angle  of  135°  with  the  positive  x-axis.  These  are  electro- 
static waves  which  can  propagate  only  for  x > 0. 


For  large  kx , e is  independent  of  kx  and  G is  proportional  to  kx.  Thus  the  integrand  is 
evanescent  as  |fcx|  — - and  Im(kx  ) < 0. 

Two  possible  contours  for  the  evaluation  of  Eq.  (20)  are  shown  in  Fig.  2 and  are  labeled 
L]  and  L2.  Both  contours  terminate  at  |/rx|  = 00  in  the  sector  of  the  fr^-plane  for  which 
0 > arg(kx)  > - 7t/3.  The  contour  £,  loops  around  the  branch  point  at  + k0  and  the  con- 
tour i2  loops  around  the  branch  point  at  — k0.  For  x > 0 both  contours  can  be  deformed  so 
as  to  pass  through  the  saddle  point  at  + ks  along  the  path  of  steepest  descent.  Thus,  as 
|x|  — °°  the  solution  obtained  from  Lx  which  we  call  £]  (x),  will  contain  a contribution  from 
the  branch  point  at  k0  (an  extraordinary  wave  with  vgx  = dw/dkx  > 0)  and  a contribution 
from  the  saddle  point  (an  electrostatic  wave  with  v^  > 0).  The  solution  obtained  from  L2, 
which  we  call  E2 , will  contain  a contribution  from  the  branch  point  at  — k0  (an  extraordinary 
wave  with  v^  < 0)  and  a contribution  from  the  saddle  point  (an  electrostatic  wave  with 
V > °>- 

For  x < 0 the  two  contours  are  deformed  as  shown  in  Fig.  3.  Contour  L2  passes 
through  the  saddle  point  now  located  on  the  imaginary  kx  axis.  However,  the  contribution 
from  this  saddle  point  is  evanescent  as  x — — 00  and  need  not  be  considered.  Thus,  only 
contributions  from  the  branch  points  need  be  considered  for  large  negative  x.  £,  (x)  will  con- 
tain a contribution  from  k0  (and  extraordinary  wave  with  vgx  > 0)  and  E2  (x)  will  contain 
contributions  from  both  branch  points  (extraordinary  waves  with  \gx  < 0).  Therefore  the 
solution  £|  (x)  applies  for  the  case  in  which  an  extraordinary  wave  is  incident  from  the  high 
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magnetic  field  side  Or  — — °°)  of  the  resonance.  The  solution  £2  Or)  contains  two  waves  in- 
cident on  the  resonance  and  in  itself  does  not  apply  to  a case  of  physical  interest.  However  an 
appropriate  linear  combination  of  £,  and  £2  can  be  formed  which  applies  to  the  case  of  an  ex- 
traordinary wave  incident  from  the  low  field  side  of  the  resonance. 
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The  asymptotic  values  of  Ei  (x)  and  £2  (x)  are  seen  to  be  linear  combinations  of  the 
various  WKB  solutions  of  the  local  dispersion  relation.  Thus,  Eq.’s  (23)  and  (24)  provide  con- 
nection formulas  linking  the  WKB  solutions  on  either  side  of  a resonant  surface.  To  determine 
the  energy  transmission  and  reflection  coefficients  we  must  determine  the  energy  flux  of  each 
of  the  WKB  solutions  in  the  asymptotic  expansions.  The  energy  flux  of  a wave  in  a weakly  in- 
homogeneous plasma  is  given  by  Stix5, 


E, 


(25) 


where  D - — - k^k*  - I +1.  In  our  case  Eq.  (25)  reduces  to 

OJ  2 


*(kx)  + -j  GUJ 


(26) 


where  k — k for  an  extraordinary  wave  and  kx  — ks  (x)  for  an  electrostatic  wave. 


A.  Wave  Incident  from  the  High  Field  Side 


With  the  definition  of  the  wave  energy  flux  given  by  Eq.  (26)  we  can  calculate  the  ener- 
gy transmission  and  reflection  coefficients  for  the  solution  £,  (x),  which  applies  in  the  situa- 
tion of  an  extraordinary  wave  incident  on  the  resonance  from  the  high  field  side.  The  ratio  of 
transmitted  to  incident  power  is  found  from  Eq.  (23a)  and  Eq.  (24a)  to  be. 


Tj  = exp  ( — 2wt)0  ). 

(27) 

This  is  precisely  the  result  predicted  by  WKB  theory.  The  ratio  of  mode  converted  electrostatic 
wave  energy  to  incident  wave  energy  is  found,  after  some  algebra,  to  be, 

FX  - (1  - 7-.  >• 

(28) 

Thus,  all  the  energy  that  WKB  theory  predicts  to  be  absorbed  is  in  fact  mode  converted  to  the 
electrostatic  wave.  If  we  now  reinsert  the  effects  of  finite  k2  by  shifting  the  x-axis,  per  the 
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discussion  following  Eq.  (20),  we  find  that  Eq.’s  (27)  and  (28)  remain  valid.  However,  the  elec- 
trostatic wave  is  heavily  damped  (cf.  Eq.  (22))  so  the  net  effect  is  the  absorption  by  the  elec- 
trons of  a fraction  F,  of  the  incident  energy.  In  this  sense  WKB  theory  is  seen  to  give  the 
correct  answer  for  Lk.k0\e/w  > 1. 


B.  Wave  Incident  from  the  Low  Field  Side 


In  order  to  obtain  a solution  which  contains  a single  incident  wave  (incident  from 
x = oo)  we  form  the  linear  combination, 

E3  (x)  = E2  (x)  -(1  - exp  ( -2irr)0  ) ]£,  (x). 

(29) 

From  Eq.’s  (24)  it  is  seen  that  E3  (x)  contains  only  a transmitted  extraordinary  wave  for 
x < 0.  From  Eq.’s  (23)  it  is  seen  that  E3  (x)  contains  an  incident  and  a reflected  extraordi- 
nary wave  and  a backscattered  electrostatic  wave  for  x > 0.  We  define  the  coefficients 
T2,  R2,  and  F2  to  be  the  fractions  of  energy  transmitted,  reflected,  and  mode  converted  to  an 
electrostatic  wave  respectively.  By  application  of  Eq.  (26)  these  coefficients  are  found  to  be 

T2  = exp  ( ~2nr)n  ). 

(29a) 

R2  = (1  - T2  )2. 


F2  = T2  (1  - T2 ) 


which  apply  for  k.  = 0.  It  can  be  verified  T2  + R2  + F2  =1. 


Inclusion  of  finite  k:  modifies  only  R2  It  is  found,  by  shifting  the  x-axis,  that  R2  is  re- 
duced by  a factor  exp  ( -2kaLk.\e/c)  for  small  but  finite  k:.  Thus,  we  see  that  for  k.  -=  0 all 
the  incident  energy  remains  in  waves,  but  if  k.K0L\e/c  — 1 a substantial  portion  of  this  ener- 
gy can  be  transferred  to  particles  by  cyclotron  damping 


•w 


VI.  ORDINARY  WAVE  AT  THE  FIRST  HARMONIC 

Our  analysis  in  this  section  parallels  that  of  the  previous  section.  Here  we  consider  the 
propagation  of  the  ordinary  mode  in  the  vicinity  of  the  first  cyclotron  harmonic.  With  k.  = 0 
a single  integral  equation  describes  the  propagation  of  the  ordinary  mode. 


Sdkx 


2n 


r2  . 2 
KX  C D 

— 1 + — 

2 2 

<0  (t) 


E,  — iL  J dk'  7 E,  (*') 


krk'\} 


(It  — oo 


40; 


exp(/A  x)  = 0, 


(30) 


where  we  have  assumed  v2  < < c2  in  evaluating  the  matrix  element  Ip.  Equation  (30)  can 
be  cast  in  the  form  of  Eq.  (18)  and  hence  yields  the  solution, 


EAx) 


r kxdkx 

4 


exp 


ikxx  + iL  J dk'G(k')lt(k') 


t ( kr ) 


(31) 


where 


and 


G(k')  = ja2v2*'2/02, 


t(kx)  — )c2c2/oj2  —1  +a2. 


Two  contours  for  the  asymptotic  evaluation  of  Eq.  (31)  are  shown  in  Fig.  4.  These  con- 
tours are  similar  to  those  of  Fig.  2 for  the  evaluation  of  Eq.  (20).  However,  in  the  present  case 
there  are  no  saddle  points,  and  the  evanescence  of  the  integrand  for  large  | A^.  | is  determined 

*x 

by  the  ikxx  term  rather  than  the  J dk'G/t  term.  Thus,  for  x > 0 we  must  end  the  contours  in 

the  upper  half  of  the  fc^-plane  and  for  x < 0 the  contours  and  branch  cuts  must  be  rotated  to 
the  lower  half  of  the  fcx-plane. 


Figure  5 shows  the  contours  for  x < 0 that  are  obtained  when  the  branch  cuts  and  con- 
tours are  rotated  in  the  clockwise  sense.  The  proper  sense  of  rotation  is  determined  from  the 
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analytic  properties  of  E.  (x).  We  see  from  Eq.  (14)  that  the  causality  condition  Im  (<u)  > 0 is 
equivalent  to  requiring  the  conductivity  elements  be  analytic  in  the  upper  half  of  the  complex 
jr-plane.  Thus,  we  seek  a solution  £.  (x)  with  this  same  property.  That  is,  the  analytic  con- 
tinuation of  E:  (x)  for  x < 0 is  obtained  from  the  solution  for  x > 0 by  increasing  the  argu- 
ment of  x from  0 to  n.  Hence,  the  branch  cuts  and  contours  of  Fig.  4 are  rotated  in  the  clock- 
wise sense  to  obtain  the  analytic  continuation  of  E . (x). 


For  large  x only  the  portions  of  the  contours  near  the  branch  points  contribute  to  the  in- 
tegrals. Thus,  we  may  write  down  the  asymptotic  expressions  for  E1;  (x)  and  £2.  (x)  by  anal- 
ogy to  Eq.’s  (23)  and  (24), 


E|  (x)  ~ k0  [1  - exp  ( -2rrT)0 ) ) r (/Tj0)exp  (7tt)0/2 )i/< +0  / 2 n-—f(k0) 


E2  (x)  — k0  [1  - exp  (2itt)0  ) ) T ( —ir)0  )exp  ( -37tt)0/2  )t //  _0  / 27t  — — e (k0  ) 


for  x > 0,  and 


E\  — [1  “ exp  ( -2ttt)0  ) ] F (/t)0  )exp  (3itt)0/2  )i l>  +0  / 2t t — — « (k0 ) 


E2  — k0  [1  - exp  (27tt/0  ) ] T ( -ir)0  )exp  ( — 5ttt)0/2  )i h _0  / 2t t e (*:„  ) 


for  x < 0,  where  k^c^/w2  =1  — 


The  determination  of  the  transmission,  and  reflection  coefficients  proceeds  as  before.  We 
find  for  the  case  of  an  ordinary  wave  incident  from  the  high  magnetic  field  side  a fraction 
T - exp  ( — 27TT)a  ),  where  t)0  = LG(kn)/  (df/dk0)  is  transmitted,  and  no  energy  is  reflected. 
This  is  equivalent  to  the  WKB  result,  Eq.  (5). 


In  the  case  of  a wave  incident  from  the  low  field  side  a fraction  of  energy 
T = exp  ( —2nr)a  ) is  transmitted,  and  a fraction  R — (1  — D 2 is  reflected. 
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In  both  of  these  cases  the  energy  contained  in  the  scattered  waves  is  less  than  the  energy 
contained  in  the  incident  wave  indicating  the  absorption  of  energy  by  the  electrons.  That  this 
absorbed  energy  goes  into  increasing  the  motion  of  the  electrons  perpendicular  to  the  applied 
magnetic  field  rather  than  the  parallel  motion  can  be  seen  by  considering  the  component  of 
canonical  momentum  parallel  to  the  magnetic  field.  For  the  case  of  purely  perpendicular  propa- 
gation, which  we  consider  here,  the  electromagnetic  wave  field  can  be  generated  by  the  vector 
potential  A = *:A  (x,  t ).  Thus  the  z component  of  canonical  angular  momentum  is  conserved, 
p.  “ mv_.  - eA  (x,  i)/c.  A(x,  t)  is  an  oscillating,  bound  function.  Thus,  v.  must  be  oscillating 
and  bound  as  well  indicating  that  on  the  average  the  electrons  do  not  gain  parallel  energy.  En- 
ergy is  fed  at  the  perpendicular  motion  by  the  e\:By/c  force  where  Bv  is  the  wave  magnetic 
field. 

If  we  allow  k.  to  be  sufficiently  large  ( k:\e/w  ^ ( kLL ) _l ) then  the  width  of  the  reso- 
nance becomes  large  enough  that  WKB  theory  applies.  In  this  case  the  transmission 
coefficients  remain  the  same,  but  the  reflection  that  occurs  for  a wave  incident  from  the  low 
field  side  with  k.  = 0 should  disappear.  Thus  the  ordinary  wave  can  be  effectively  absorbed  at 
the  first  cyclotron  harmonic  (w  = il  (x))  regardless  of  whether  it  is  incident  from  the  high  or 
low  field  side  of  the  resonance. 

VII.  ACCESSIBILITY  OF  CYCLOTRON  RESONANCES  AND  APPLICATIONS  TO  HEATING  IN 
TOROIDAL  DEVICES 

In  the  preceeding  sections  of  this  paper  we  found  two  situations  in  which  electromagnetic 
wave  energy  could  be  absorbed  directly  at  a cyclotron  harmonic.  These  situations  were:  an  or- 
dinary wave  incident  on  the  first  cyclotron  harmonic,  and  an  extraordinary  wave  incident  on 
the  second  harmonic.  It  was  also  found  that  there  could  be  appreciable  energy  absorption  re- 
gardless of  whether  the  wave  was  incident  on  the  resonance  from  the  high  or  low  magnetic 
field  side 
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The  accessibility  conditions  of  the  first  cyclotron  harmonic  for  the  ordinary  wave  are 
quite  simple.  If  k.  =0  all  that  is  required  is  that  the  plasma  density  be  low  enough  so  that 
(Ope  < co2  = ( l2  (c.f.  (Eq.  (1)).  It  is  anticipated  that  oj2(,  < llj  for  most  tokamaks  so  the 
ordinary  wave  shows  great  promise  for  plasma  heating  purposes.  Allowing  for  nonzero  k.  the 
accessibility  condition  becomes  (t»2e/oj2)  < the  minimum  of  1 and  2 (1  — k}c2/uj2). 

In  order  for  the  second  harmonic  to  be  accessible  to  the  extraordinary  wave  we  must 
have  w2 /ia 2 = yco2/n2  < y (c.f.  Eq.  (2)),  or  ca2  < 2 ii  2.  Thus  the  second  harmonic  is 

available  for  heating  in  higher  density  plasmas.  Inclusion  of  finite  k:  alters  this  condition 
slightly,  ca 2 < 2 SI  2 ( 1 —k2c2/ ca2). 

In  Fig.  6 the  three  proposed  methods  of  heating  using  waves  in  the  electron  cyclotron 
frequency  range  are  summarized.  The  ray  paths  through  the  parameter  plane 
(ca2(,  Or)/ca 2,  fl2(.v)/ca2)  are  traced  for  the  ordinary  wave  incident  on  the  first  harmonic,  the 
extraordinary  wave  incident  on  the  second  harmonic,  and  the  extraordinary  wave  incident  on 
the  upper  hybrid  resonance.  The  paths  of  incidence  from  both  high  and  low  magnetic  field 
sides  are  shown  for  the  two  cases  of  direct  absorption  at  a cyclotron  resonance  considered  in 
this  paper. 
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Appendix 

ASYMPTOTIC  EVALUATION  OF  E+  (*) 

The  contours  for  the  evaluation  of  Eq.  (20)  are  shown  in  Fig.  2.  Branch  cuts  in  the 
Av-plane  originate  from  each  of  the  solutions  of  t (A)  =0,  ( = ± k0).  These  branch  cuts  are 
associated  with  different  possible  paths  of  integration  for  determining  the  function. 


M(kx)  = f dk'G(k')/e  (A  ). 


(Al) 


A.  Saddle  Point  Contribution 


The  contour  L|  loops  around  the  branch  point  at  kg  and  for  x > 0 both  the  contour  and 
branch  cut  are  made  to  pass  through  the  saddle  point  along  the  path  of  steepest  descent.  The 
contour  passes  through  the  saddle  point  twice.  The  value  of  the  integrand  on  the  upper  (lower) 
leg  of  the  contour  is  determined  by  choosing  the  path  of  integration  in  M(kx)  to  pass  above 
(below)  the  singularity  at  e ( k0  ) = 0. 

To  evaluate  the  saddle  point  contribution  we  define  F to  be  the  exponent  of  the  in- 
tegrand. 


F(kx ) = ikxx  + iLM(kx), 


(A2) 


which  we  will  Taylor  expand  to  second  order  in  kx  at  the  saddle  point,  kx  = Avs.  The  contri- 
bution from  the  saddle  point  is  therefore  given  approximately, 


Es  — € 1 (Avv ) exp  (F(kxs ) ) f rfA'exp 

2 7T  U 


\F"(kxs)(k  — kxs ) 2 


where  the  vanishing  of  the  first  derivative  of  F(ky  ) defines  the  saddle  point. 


dF 

~ = ix  + iLG(kxs  )/t  (kxs ) = 0. 


(A3) 


Equation  (A3)  can  be  rearranged  in  such  a way  so  that  it  is  seen  to  be  identical  to  the  local 
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dispersion  relation  We  define  the  solution  of  Eq.  (A3)  to  be  = kxs  ( v ) . 


r 


* 


The  value  of  F(kxh  (x) ) must  be  determined  for  both  legs  of  the  contour, 

(kxs(x))  = ik  vv  (x)x  + iLP  f dk'G(k')/t  (k‘)  ± ttVo. 

*0 

where  r)0  = LG(ka)  [de  (ko)/Bk0  ] and  4-  ( — ) applies  to  the  upper  (lower)  leg  of  the  con- 
tour. By  simple  manipulation,  and  use  of  Eq.  (A3)  we  find, 

X 

F+  (kxs  (x))  — Pi  f dx'kxs(x')  ± 7r  r/  () . 

*0 

Thus,  the  saddle  point  contribution  will  match  on  to  the  WKB  solution  of  an  electrostatic  wave 
obtained  from  the  local  dispersion  relation. 

Upon  expanding  F(kx)  to  second  order  in  kx  — A(s,  performing  the  indicated  integra- 
tions, and  summing  the  contributions  from  both  legs  of  the  contour,  we  obtain  the  saddle 
point  contribution  to  F.\  (.v)  for  x > 0, 

F.c  | = — sinh  w»)0exp  [/ («  -4-  3 7t/4  ) 1 x (2/  ( -nr  e 0 («.v  4-  Z.G)/dAu  ) ) 1 2 

(A  5) 

where  all  quantities  are  defined  in  the  text.  T he  saddle  point  contribution  to  F2  ( v)  is  found  in 
a similar  manner  to  be  — Et]  (a  ). 

B.  Branch  Point  Contribution 

Let  us  define  // ( A v ) in  the  following  way, 

*«  *, 

iL  J dk'G/t  = i/l(kx  ) 4 n,n  f dk’lk  — k ) -1, 

()  0 

where. 


I, 
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iH(kr)  = iL  f dk 
Jo 


G(k') 


( k ') 


G(k„ ) 

t'(k0)(k'  - k0 ) 


is  an  analytic  function  of  kx  near  the  branch  point  k„.  We  now  transform  variables  of  integra- 
tion by  letting, 


k — k0  = z/exp  ( —3  in/ 2 )/x. 

(A6) 

The  contour  of  integration  in  the  u- plane  is  shown  in  Fig.  7.  The  transformation,  Eq.  (A6),  is 
such  that  2n  > argiu)  > 0.  Thus  for  x > 0 and  |x|  — w e obtain, 


F.m  = — k0z\v\ik0x  + z/y,  (*0)  - (T)0ln  (kax) 


+ 7rr)0/2!c  (r)0 )/  (dt/dk„), 


(A7) 


where 


c(v0)  = / duii  'expl-w  + iri0lnu]. 


c(v0)  = n — exp  ( — 2ny)0  ) ) r (/T)0  ). 
For  x < 0 it  is  necesary  to  let  x = | jc|  p ,7T. 


In  a similar  manner  the  branch  point  contributions  for  £2  can  be  obtained.  However, 
there  are  contributions  from  the  portions  of  the  contour  L2  that  pass  near  the  branch  points  at 
kn.  These  contributions  are  found  to  be  linearly  proportional  to  the  asymptotic  value  of  E,  (x). 
The  final  asymptotic  formulas  are  given  in  Eq.  (23)  and  Eq.  (24). 
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Figure  1 — Allis  diagram  for  electromagnetic  wave  propagation 
with  fixed  kz  ( kz  II  B)  in  a cold  plasma.  The  solid  line  is  the 
upper  hybrid  resonance.  The  dashed  lines  are  various  cut-offs. 
In  regions  containing  the  letter  K(O)  the  extraordinary  (ordi- 
nary) wave  propagations. 


2 


Figure  3 — The  deformation  of  the  two  contours  of 
Fig.  2 for  the  evaluation  of  Eq.  (20)  for  x < 0. 


